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Collective quadrupole and octupole states are described in a series of Sm and Gd isotopes within 
the framework of the interacting boson model (IBM), whose Hamiltonian parameters are deduced 
from mean field calculations with the Gogny energy density functional. The link between both 
frameworks is the (/?2/93) potential energy surface computed within the Hartree-Fock-Bogoliubov 
framework in the case of the Gogny force. The diagonalization of the IBM Hamiltonian provides 
excitation energies and transition strengths of an assorted set of states including both positive 
and negative parity states. The resultant spectroscopic properties are compared with the available 
experimental data and also with the results of the configuration mixing calculations with the Gogny 
force within the generator coordinate method (GCM). The structure of excited 0’*' states and its 
connection with double octupole phonons is also addressed. The model is shown to describe the 
empirical trend of the low-energy quadrupole and octupole collective structure fairly well, and turns 
out to be consistent with GGM results obtained with the Gogny force. 

PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60.Jz 


I. INTRODUCTION 

The study of the equilibrium shapes and the corre¬ 
sponding excitation spectra of atomic nuclei is one of 
the recurrent themes in nuclear structure physics. Most 
of the deformed medium-heavy and heavy nuclei exhibit 
reflection-symmetric ground states. However, in some re¬ 
gions of the nuclear chart, there is an onset of reflection- 
asymmetric shapes driven by specific shell effects. In 
quadrupole deformed nuclei, a characteristic feature of 
octupole deformation is the alternating-parity rotational 
band formed by the even-spin positive parity states and 
alternating odd-spin negative-parity states, connected 
with each other by enhanced electric dipole transitions 

ra¬ 
in the framework of the spherical shell model, octupo- 
larity arises as a result of the coupling between the {l,j) 
orbitals in a major shell and the unique-parity (1-1-3, j-|-3) 
intruders from the next major shell. Within this context, 
illustrative examples are the rare-earth nuclei with the 
proton number Z « 56 and the neutron number N ks S8 
as well as the light actinides with Z « 88 and N « 134. 
In the light actinides case, the coupling of both neutron 
(i-e-, 1^9/2 and Ojis/a) and proton (i.e., I/7/2 and 0*13/3) 
single-particle states leads to octupole deformed ground 
states [lira- A recent Coulomb excitation study has re¬ 
vealed, for the first time, unambiguous evidences of static 
octupole deformation in ^^"‘Ra [3]. 

In this work, we study the impact of octupole corre¬ 
lations on the ground state and the associated low-lying 
collective spectra of the nuclei and 

We consider both quadrupole and octupole degrees of 
freedom. The selected nuclei belong to a region of the 
nuclear chart where octupole correlations are expected to 
play an important role and therefore, represent a valuable 


testing ground for the considered theoretical approxima¬ 
tions. Indeed, the experimental observation of octupole 
correlations at medium spin, as well as the crossing of the 
octupole and the ground-state bands, point to the coexis¬ 
tence of reflection symmetric and asymmetric structures 
in both ^^°Sm [3] and ^^®Sm 0- From the experimen¬ 
tal point of view, four low-lying negative-parity bands 
have already been identified in ^^^Sm [6]. The emerging 
pattern of excitations, suggests a complex shape coexis¬ 
tence in this nucleus. Moreover, the nucleus ^^^Sm has 
been identified [7] as an example of the X(5) critical point 
symmetry [8]. The nature of many low-lying excited O’*" 
states in rare-earth nuclei has also attracted much at¬ 
tention. For example, thirteen excited O’*" states have 
already been identified for ^®®Gd [3]. Within the spdf- 
IBM framework, many of the observed O’*" states have 
been attributed to the coupling of two octupole phonons 

m- 

Keeping in mind the experimental findings mentioned 
above, it is interesting and timely to consider a system¬ 
atic analysis of the quadrupole-octupole collectivity in 
rare-earth nuclei. The breaking of reflection symme¬ 
try and the associated low-lying negative-parity states 
have been addressed using various theoretical frame¬ 
works: self-consistent mean-field [IIHID], algebraic [HI¬ 
ES] , collective phenomenological [2M33] , and cluster [34l - 
I36j models. A large number of calculations for nuclei 
with static and/or dynamical octupole deformations have 
already been reported [12 fTilEm [371j33| . In particu¬ 
lar, the nuclear energy density functional (EDF) frame¬ 
work, both at the mean-field level and beyond, provides 
a reasonably accurate description of the properties of the 
negative- and positive-parity states all over the nuclear 
chart [40]. Both non-relativistic [4TH43] and relativistic 
[331133] EDFs have already been applied in both mean- 
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field and beyond mean-field studies of medium-heavy and 
heavy mass nuclei. The description of the excitation 
spectra and transition rates requires the inclusion of dy¬ 
namical (i.e., beyond mean-field) correlations associated 
with the restoration of the broken symmetries and/or 
fluctuations in the collective parameters (i.e., generating 
coordinates) [HI [HJ 001 05] . Within this context, the 
projection of the intrinsic (i.e., symmetry-broken) states 
onto good parity ones as well as the corresponding con¬ 
figuration mixing, in the spirit of the two-dimensional 
generator coordinate method (GCM) [37], have been con¬ 
sidered recently for nuclei in the rare-earth region using 
the quadrupole Q 20 and octupole Q 30 moments as gen¬ 
erating coordinates [19]. For recent GCM study, based 
on ( 530 -constrained mean-field states, the reader is also 
referred to Ref. [15] . 

In this work we first carry out ((520 j Q 3 o)"Constrained 
Hartree-Fock-Bogoliubov (HFB) calculations based on 
the Gogny-EDF [43]. Such calculations provide us with 
the corresponding (axially symmetric) mean-field poten¬ 
tial energy surfaces (PES). Subsequently, in order to 
obtain the spectrum and wave functions of the excited 
states, we employ the interacting boson model (IBM) 
[15] . The essence of our method is to determine the pa¬ 
rameters of an appropriate IBM Hamiltonian by calcu¬ 
lating the associated bosonic PES so that it matches the 
Gogny-HFB PES. The IBM Hamiltonian resulting from 
our fermion-to-boson mapping procedure is then used in 
spectroscopic calculations. A similar mapping has been 
used in previous studies of low-lying quadrupole states 
[501 - 152] and shape coexistence [53] . Recently, the method 
[15] has been extended to describe quadrupole-octupole 
correlations and shape transitions in the light actinide 
and rare-earth regions [HES] based on the relativistic 
DD-PCI EDF. 

The same Gogny-EDF can be used along with be¬ 
yond mean field techniques to restore the broken reflec¬ 
tion symmetry and compute the properties of the lowest 
lying negative parity state. The excitation energy and 
transition strengths, when compared with the IBM num¬ 
bers, can be used as a benchmark to test the consistency 
of the mapping procedure. Therefore, one of the goals 
of this study is to assess the fermion-to-boson mapping 
methodology in the description of spectroscopic proper¬ 
ties in rare-earth nuclei. We compare the IBM spectra 
and transition rates with previous Gogny-GGM calcu¬ 
lations for the same Sm and Gd nuclei m as well as 
with available experimental data. Here, we also refer the 
reader to the previous IBM study based on the relativis¬ 
tic mean-field (RMF) approximation [SS]. We have used 
the DIM [55] parametrization of the Gogny-EDF, which 
was originally designed to better describe nuclear masses. 
It has been shown [BETHSO] that the DIM parameter 
set essentially retains the same predictive power as the 
standard and thoroughly tested Gogny-DlS [61] one. We 
have also performed a selected set of calculations based 
on the DIS parametrization in order to examine the ro¬ 
bustness of our predictions with respect to the particular 


version of the Gogny-EDF employed. However, as the 
corresponding HFB [15] and IBM results are quite simi¬ 
lar, in the present paper we will only focus on calculations 
based on the DIM parameter set. 

The paper is organized as follows. In Sec.]!^ we briefly 
outline the HFB-to-IBM mapping procedure. Next, in 
Sec. HI we discuss the systematics of the {1320, Pso) [52] 
PESs obtained for the considered nuclei as well as the 
parameters of the IBM Hamiltonian. The results of 


the spectroscopic calculations are discussed in Sec. IV 


First, in Sec. jlV A] we present the systematics of the low- 
energy spectra and the reduced transition probabilities in 
i46-i56gj^ and ^^®“^®®Gd. We will compare with avail¬ 
able experimental data as well as with results obtained 
within the Gogny-GGM approximation [T5]. Next, in 
Sec. |IVB| we further illustrate the predictive power of 
the mapped IBM model with a detailed discussion of the 
spectroscopic properties for ^®°Sm (a soft nucleus along 
the quadrupole and octupole directions) and ^®®Gd (a 
strongly quadrupole deformed nucleus). In order to ob¬ 
tain some insight into the nature of the excited O'*" states 
in the studied nuclei, their systematics is discussed in 
Sec. jlV C| In Sec |IVD[ we discuss the IBM correlation 
energies and compare them with Gogny-GGM results. 
Finally, Sec. [V] is devoted to some concluding remarks 
and work perspectives. 


II. FRAMEWORK 

In this section we briefly outline the HFB-to-IBM map¬ 
ping scheme [55] . Our starting point is a set of axi¬ 
ally symmetric ((5207 Q 3 o)-constrained Gogny-HFB cal¬ 
culations m- They provide us with the corresponding 
mean-field potential energy surfaces (MFPESs) and the 
HFB states |4’((520) Q 30 )) for the nuclei and 

i4®“i®®Gd. For simplicity, both the quadrupole Q 20 and 
the octupole Qso moments are then translated into the 
standard P 2 and mean-field deformation parameters. 

Subsequently, the MFPESs obtained are mapped into 
their bosonic counterparts, i.e., the IBM potential energy 
surfaces (IBMPESs). This procedure allows us to deter¬ 
mine the parameters of the IBM Hamiltonian used in the 
spectroscopic calculations. The IBM Hamiltonian is con¬ 
verted into a potential energy surface by means of a set 
of coherent bosonic states and this IBM-PES is what is 
used to match the Gogny-HFB PES [55]. Note that the 
MFPESs correspond to the total HFB energies, i.e., nei¬ 
ther mass parameters nor zero point (rotational and/or 
vibrational) quantum corrections are included. 

The description of the quadrupole and octupole defor¬ 
mations as well as the positive- and negative-parity states 
within the IBM framework requires both positive- and 
negative-parity bosons. Here, one assumes that the low- 
lying positive-parity states are reasonably well described 
by the pairs of valence nucleons associated to the s and d 
bosons, respectively. On the other hand, negative-parity 
states are assumed to be described by the coupling to 
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octupole / bosons [53] . Therefore, our entire IBM model 
space comprises the s, d and / bosons. For simplicity, we 
do not distinguish between proton and neutron bosons. 
A more complete description of the low-energy collective 
states would require the inclusion of the dipole p boson 
that could be associated to the spurious center-of-mass 
motion [33] or to the giant dipole resonance [M]. This, 
however, lies out of the scope of the present paper and is 
left for future work. 

The sdf Hamiltonian used is given by 

H = e-dfi-d + e/fi/ + K 2 Q 2 ■ Q 2 + • Ld + 

where the hrst (second) term stands for the number op¬ 
erator for the d (/) bosons with (e/) being the single d 
(/) boson energy relative to the s boson one. The third 
term represents the quadrupole-quadrupole interaction 
with strength K 2 - The quadrupole operator is given as 

Q 2 = d + s + Xddld^ X +Xff[P X (2) 

where Xdd and Xff are parameters. The forth term in 
Eq. Q is the rotational one relevant for the sd space. In 
this case, the angular momentum operator Ld reads 

La = '/I0[d1’ X (3) 

The last term in Eq. 0 is the octupole-octupole inter¬ 
action with the strength parameter K 3 . The octupole 
operator takes the form 

Qs = f + f'’s + Xdfld^ X / -h/''' X (4) 

with Xdf being a parameter. 

Note, that Eq. Q does not represent the most gen¬ 
eral form for the sdf Hamiltonian. The present form 
has already been used in previous phenomenological IBM 
studies which have confirmed its suitability to describe 
the available experimental data. The Hamiltonian 
of Eq. Q can be derived from a microscopic octupole- 
octupole interaction between proton and neutron bosons 
by mapping the totally symmetric state in the IBM-2 
space onto the equivalent one in the IBM-1 space [55]. We 
neglect the dipole-dipole interaction term Ld ■ Lf (with 
Lf = '/^[df X /]^^^), because it has been shown [35] to 
be of little relevance for low-energy states. 

The IBMPES is calculated as the expectation value of 
the Hamiltonian Eq. 0 in the boson condensate state 



\4>) = with -b Ml + / 33 /d. 

\/I\B 

(5) 

where Nb{= ris + Ud + rif) and |—) denote the total 
number of bosons (i.e., half the number of valence nucle¬ 
ons [63] 1 and the inert core, respectively. In the present 
study, the doubly-magic nucleus ^^^Sn is assumed to be 
the inert core. Therefore, Nb runs from 6 to 12 (7 to 
13) in (^^®“^^®Gd). For the quadrupole case 

(A = 2) the bosonic j32 and fermionic /32 deformations 


can be related as $2 = C' 2/32 [55], with C 2 being a coef¬ 
ficient. Here, as in previous works [Ml El], we assume 
that dz = C' 3 / 33 , with Cz being an additional coefficient. 

In order to reduce the computational effort, it has 
been customary in many of the previous phenomenolog¬ 
ical IBM calculations to restrict the maximum number 
of / bosons to 71™“^ = 1 in the diagonalization of the 
IBM Hamiltonian. However, as shown in the next sec¬ 
tion, the microscopic PESs may exhibit a sizable ground 
state octupole deformation which requires a larger num¬ 
ber of / bosons in our IBM calculations. Therefore both 
positive- and negative-parity bosons are treated on an 
equal footing. As a consequence, a truncation on 
is not used and the number of / bosons can run from 0 
to Nb- This also holds true for the s and d bosons. Let 
us also mention, that previous phenomenological stud¬ 
ies (e.g., [571155]) have also suggested the need of more 
negative-parity bosons for a better description of the ex¬ 
perimental data. 

The analytic IBMPES reads 


E{P2j3) 


Nb 

1 + 4 " + M 



Nb{Nb-1) 




2 



( 6 ) 


with 

7 

e's = 5k2 - 7k 3, = Cd + 64 + (1 + - gXd/«3 

5 

and Cf = ef — -X //'«2 + (1 + Xd/)^ 3 - (J) 


The IBMPES E{( 32 ,^ 3 ) is specified by the parameters 
of the Hamiltonian in Eq. 0 plus the coefficients C 2 
and C 3 . We have determined those parameters by fit¬ 
ting the IBMPESs to the Gogny-DIM MFPESs using 
the same procedure as in Ref. m- Let us remark that, 
even though a simplihed Hamiltonian Eq. Q is consid¬ 
ered, there is still a larger number of parameters to be 
determined, as compared to the sd IBM system. There¬ 
fore, rather than trying to fit all the parameters at once, 
we first determine the ones relevant for the sd space (e^, 
1 ^ 2 , Xdd, C 2 and 4) s-iid then those associated to the / 
space as well as the ones associated with the coupling 
between the two spaces (e/, K 3 , Xff, Xdf and C 3 ). The 
Ld ■ Ld term in Eq. 0 does not contribute to the PESs, 
and therefore its strength K 2 is determined independently 
by comparing the fermionic and bosonic cranking mo¬ 
ment of inertia (see Ref. m for details). The (fermionic) 
Thouless-Valatin [69] moment of inertia for the 2^ state 
reads 


Itv = ^/E-y. 


( 8 ) 
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where stands for the 2 ^ excitation energy obtained 
from the self-consistent cranking calculation with the 
constraint {J^) = \/J{J + 1), where represents the 
X component of the angular momentum operator. On 
the other hand, the IBM moment of inertia is computed 
using the coherent state |^(/ 3 , 7 )) and the Schaaser-Brink 
m expression 


^Itrm = lim 


1 {4'{P,l)\Lx\(l){/3,'y)) 

( 0 (^, 7 )|(;()(/ 3 , 7 )) 


(9) 


with oj being the cranking frequency. 

Having the parameters e^(= Cd — 6 K 2 ); ^ 2 , Xdd and C 2 
already determined from the fit of the IBMPES to the 
MFPES in the sd space, the IBM moment of inertia in 
Eq. @ depends only in the parameter K 2 whose value is 
determined so that Iibm is equal to the Itv value at the 
energy minimum. 

From the diagonalization of the S(i/-IBM Hamilto¬ 
nian, we have obtained both the energies and wave func¬ 
tions of the spectrum which are labeled by total spin 
and parity quantum numbers. We have used the com¬ 
puter program OCTUPOLE [7T]. The reduced elec¬ 
tromagnetic transition probabilities B{E\\J —>■ J') = 
|(j/||^(£;a)||j)| 2 /( 2 j (-y _ 1^2,3) are then com¬ 

puted using the resulting IBM wave functions. Here, J 
(J') denotes the spin for the initial (final) state. Of par¬ 
ticular interest for the present study are the dipole El, 
quadrupole E2, and octupole E3 transition probabilities 
defined in terms of the operators 


= ei[d) X f + d'f X /](!) 

(10) 

^ g2Q2 

(11) 


(12) 


where Q 2 and Q 3 are the quadrupole and octupole op¬ 
erators appearing in the IBM Hamiltonian and e^’s are 
boson effective charges which are kept constant for all 
the considered nuclei. Their values are taken from previ¬ 
ous phenomenological IBM studies (ei = 0.01 eb^/^ [HH] . 
62 = 0.13 eb [ 68 ] and 63 = 0.099 eb^/^ |23|). It has been 
shown that they provide a reasonable overall description 
of the experimental data. However, they are not the ones 
derived microscopically. Therefore, in the following dis¬ 
cussions, one should always keep in mind that there is 
some extra freedom in the overall scale of the calculated 
IBM transitions. 


III. MEAN-FIELD POTENTIAL ENERGY 
SURFACES AND THE PARAMETERS OF THE 
IBM HAMILTONIAN 


In this section, we discuss the systematics of the MF- 
PESs and IBMPESs as well as the parameters of the IBM 
Hamiltonian obtained along the lines described in Sec. [TT| 
The axially symmetric Gogny-DIM MFPESs are 
shown in Fig. 0 for and The MF¬ 

PESs of some of the Sm isotopes have already been pre¬ 
sented in Ref. m as illustrative examples. However, for 


the sake of completeness, in the figure we have included 
all the MFPESs both for Sm and Gd nuclei. For the 
sake of presentation, the plots in the figure correspond 
to —0.3 ^ P 2 ^ 0.5 and 0.0 < < 0.2 as well as to 

an energy range of 5 MeV from the absolute minimum. 
We have tested, that the previous ranges are enough to 
describe the considered low-energy collective states and 
used them to build our IBM Hamiltonian. 

A spherical reflection-symmetric ground state is pre¬ 
dicted for the nuclei ^^®Sm [panel (a)] and ^"^®Gd [panel 
(g)], respectively. On the other hand, the MFPESs be¬ 
come soft for isotopes with neutron numbers N = 86 and 
88 , indicating that the Gogny-HFB approximation can 
only be considered as a valuable starting point in such 
nuclei but beyond mean-field correlations should be taken 
into account [19]. Moreover, the N = 88 isotopes exhibit 
the softest MFPESs with a shallow minimum at a non¬ 
zero jds value. One also sees that the MFPESs become 
steeper along the (is direction for isotopes with N ^ 90. 
Similar trends have been found up to A = 88 in previous 
RMF calculations [371 [^, based on the EDFs PKl [73] 
and DD-PCl [73], respectively. However in those calcu¬ 
lations, the octupole minima are more pronounced than 
ours. In fact, the previous study with the relativistic 
functional DD-PGl [35] suggested that the potential en¬ 
ergy surface is much more softer along /?3 direction. The 
same trend was found for isotopes with N ^ 90. 

As already discussed in Ref. m, there is no essential 
difference between the overall topology of the MFPESs 
obtained with the Gogny-DIM and Gogny-DlS EDFs. 
However, at a quantitative level, the latter provides MF¬ 
PESs with slightly deeper absolute minima than the for¬ 
mer. Nevertheless, such a difference turns out to be too 
small to significantly affect neither the IBM parameters 
nor the energies and wave functions of the excited states. 
With this in mind, in what follows only results based on 
the Gogny-DIM EDF will be discussed. 

In Fig. we have depicted the (mapped) IBMPESs. 
First, we observe that they are much flatter than the 
HFB MFPESs (see. Fig. [^. This is a common feature 
of the IBM framework already found in previous stud¬ 
ies [li [5D]. The reason is that IBM’s model space is 
rather limited and only comprises pairs of valence nucle¬ 
ons. This leads to flat IBMPESs for larger deformations. 
However, one should keep in mind that within the con¬ 
sidered fermion-to-boson mapping, the topology far away 
from the absolute minimum is not relevant as long as 
we restrict our analysis to the low-lying collective states. 
Hence, we only focus on reproducing the curvatures of 
the Gogny-DIM MFPESs in the neighborhood (a 5 MeV 
window) of the absolute minimum, along both the j32 and 
directions. 

Second, we note that, for N = 86 and 88 isotopes, the 
MFPES predicts a shallow absolute minimum at non-zero 
/33 values [Fig.[^ while in the corresponding IBMPES the 
absolute minimum is found at /33 = 0 [Fig. [^. However, 
as the depth of this absolute minimum in the MFPESs 
differs by at most tens of keVs from the saddle point 
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FIG. 1. (Color online) Axially symmetric (/32, Ps) potential energy surfaces for the nuclei and calculated 

within the constrained Gogny-HFB approach based on the DIM parametrization. The contour lines join points with the same 
energy (in MeV) and the color scale varies in steps of 100 keV. The energy difference between neighboring contours is 0.5 MeV. 
These ((32, Pz) energy surfaces are symmetric with respect to the Ps = 0 axis. Thus, they are only plotted for ps ^ 0. For each 
nucleus the absolute minimum is identified by an open circle. 


on the ^3 = 0 axis, we assume that the discrepancy of 
the absolute minimum point, that is not deep enough in 
energy, between the MFPES and the IBMPES is not of 
crucial importance for the final result. 

Bearing those in mind, the IBMPESs in Fig. [^closely 
follow, for each of the considered nuclei, the basic topol¬ 
ogy as well as the overall systematic trend of the Gogny- 
HFB ones shown in Fig. [T] 

In Fig. the IBM parameters for the considered Sm 
and Gd nuclei are plotted as functions of neutron num¬ 
ber. As can be observed in panels (a) and (b), the single 
d (cd) and / (e/) boson energies decrease as functions of 
neutron number. From a microscopic point of view, as al¬ 
ready discussed in the context of the sd IBM-2 [nnnES] 
model, the decrease of €d could be related to the coupling 
of the unperturbed d boson with other types of bosons 
not yet explicitly included in the model space. Alterna¬ 
tively, when one derives the form of the IBM Hamiltonian 
in Eq. ([^ from a general sdf IBM Hamiltonian, several 


two-body terms of the general IBM Hamiltonian, that are 
reduced to the kinetic energies of d and / bosons multi¬ 
plied with the boson-number dependent factors, are ab¬ 
sorbed in Cd and Cf, thereby making the parameters vary 
significantly with boson number m- 

The coupling strength of the quadrupole-quadrupole 
interaction K 2 , shown in panel (c), is almost constant. 
A similar trend has been found in the IBM study based 
on the RME approximation . A sudden change is ob¬ 
served in the parameter Xffj plotted in panel (d), around 
N = 88 and is correlated with the significant change ob¬ 
served in the MFPESs (see. Fig. [^. On the other hand, 
at variance with our previous sd IBM study in the same 
mass region m, the parameter Xdd [panel (e)] is rather 
constant. Gompared to the quadrupole-quadrupole cou¬ 
pling K 2 [panel (c)], the strength of the octupole-octupole 
interaction [panel (f)] exhibits a gradual decrease with 
increasing neutron number. 

In panel (g) of the same figure, we have plotted the 
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FIG. 2. (Color online) The same as Fig. but for the mapped IBM potential energy surfaces. 


strength K 2 of the Ld ■ Ld term Eq. ([^. Its negative 
value, for all the studied nuclei, leads to the lowering of 
the positive-parity yrast states m- Note that is not 
considered for the spherical nuclei and ^^®Gd. As 

shown below, the experimental spectra for these nuclei 
do not exhibit a rotational-like structure and, therefore, 
there is no obvious reason for introducing the Ld ■ Ld 
term in the corresponding calculations. The parameters 
Xdf [panel (h)] exhibits a pronounced isotopic dependence 
with a maximum around N = 88 — 90 which correlates 
well with the octupole softness of the MFPESs around 
the same neutron numbers. Both the C 2 [panel (i)] and 
Ca [panel (j)] coefficients change smoothly with neutron 
number [55]. 


IV. SPECTROSCOPIC CALCULATIONS 


In this section, we discuss the results of the calcu¬ 
lations with the IBM Hamiltonian for and 

i4®“i®®Gd. First, in Sec. IV A the systematics of the 


low-energy spectra and the reduced transition probabil¬ 
ities in 446-i56gj^ i48-158Qj^ jg addressed. Next, in 


Sec. IV B[ the spectroscopic properties predicted for the 
nuclei and ^®®Gd are discussed in detail. The sys- 


tematics of the excited 0+ states is presented in Sec. IV C 


Finally, in Sec |IVD| ground state correlation energies are 
discussed. 


A. Systematics of the low-energy spectra and the 
reduced transition probabilities in and 

148-158Gd 

In Figs. [^ and [^ the low-energy positive- and negative- 
parity yrast states, as calculated with the mapped sdf 
IBM Hamiltonian are plotted for the nuclei 446-i56gj^ 
and 448 -i58QjJ theoretical results are compared 

with the available experimental data taken from the 
NNDC compilation Since our predictions for Sm 

[panels (a) and (b)] and Gd [panels (c) and (d)] isotopes 
are rather similar, we mainly discuss the former. 

The lowering of the energies with increasing neutron 
number N is consistent with a shape transition (see, 
Fig.[^ to a strongly quadrupole deformed configurations. 
Indeed, the ratios Ri /2 = E{L^)/E{2'^)=2.'i2> and 2.38 
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FIG. 3. (Color online) The parameters of the sdf IBM Hamiltonian H in Eq. Q, as well as the proportionality coefficients C 2 
and Ca, are plotted as functions of the neutron number for the considered nuclei. The parameters Xdd, X//i Xdf, C 2 and C 3 
are are dimensionless. 


obtained for are both close to the vibrational 

limit while the theoretical (experimental) R 4/2 values for 
the transitional nuclei are 2.82 (2.31) and 2.91 

(3.01), respectively. Our calculations predict a more pro¬ 
nounced rotational character for ^^°Sm than expected 
from the experiment. On the other hand, it is remark¬ 
able that the R 4/2 value for the ^®^Sm is exactly the same 
as the X(5) one [5]. For the heavier isotopes, our IBM 
calculations predict well developed rotational bands. For 
example, in the case of we have obtained the 

ratios i?4/2=3.21 and 3.25, respectively. The theoreti¬ 
cal results agree reasonably well with the experimental 
ones except for the lightest isotopes where the energies 
of the higher spin states are overestimated. The rea¬ 
son for the overestimation could be the too restricted 
model space and/or Hamiltonian of the IBM that is not 
rich enough as to reproduce the peculiar topology of the 
Gogny-EDF MFPES for the lightest isotopes. We recall 
that the term is not included in and ^^®Gd 

as it is of little importance for these spherical nuclei m- 
One could introduce this term phenomenologically to fix 
the overestimation, which is however out of scope of the 
present work. 

The J'" = I“,..., 9“ states, plotted in Eig. display 


features characteristic of the octupole collectivity. Ex¬ 
ception made of the 3“ states, their excitation energies 
decrease sharply for SA < N < 90. At variance with 
the experimental data, the theoretical excitation ener¬ 
gies increase for N > 90 which correlates well with the 
diminishing of the octupole minimum depth observed in 
the MFPESs (see, Eig. [^. In both isotopic chains, the 
3/ state is lower in energy than the 1/ one. We have also 
found a near degeneracy for the 1/ and the 5/ states for 
N ^ 88 — 90. This octupole vibrational feature becomes 
more apparent for the lighter isotopes. 

In Eig. 1^ we have compared the excitation energies 
of the lowest 1/ states with the ones obtained in the 
framework of a two-dimensional GGM calculations [TS] 
also with the Gogny-DIM EDE. The predicted IBM and 
GGM values are quite similar for 84 < < 88. In the 

case of the Sm isotopes both the GGM and IBM exci¬ 
tation energies increase with increasing neutron number 
though the former exhibit a more pronounced change 
than the latter. Similar results are obtained for Gd iso¬ 
topes, exception made of the fact that the smallest 1“ 
excitation energy is found at = 88 {N = 90) in the 
GGM (IBM) calculations. 
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Neutron Number Neutron Number 

FIG. 4. (Color online) The energy spectra of the lowest-lying even-spin positive-parity states up to = 10^ for the considered 
Sm and Gd isotopes. All the experimental data are taken from the NNDC compilation m- 


We have studied the quantity 

^ [E{J+l)-E{J)]-[E{J)-E{J-l)] 

^ ’ E{2+) ’ 

which is sensitive to the splitting between the positive- 
and negative-parity members of a rotational band. In 
Eq. ( |I^ , E{J) stands for the excitation energy of the 
J = 0^, 1“, 2+, 3“,... state. Note that, for an ideal 
alternating-parity band, we would obtain an equal en¬ 
ergy splitting between the positive- and negative-parity 
states differing by AJ = 1. This, in turn, would lead to 
S{J) « 0. On the other hand, a non-zero S{J) value in¬ 
dicates a deviation from a pure alternating-parity band. 

In Fig. [^we have plotted as a function of the 

spin J, for 150452,i56gj^ which are taken as representa¬ 
tive examples. The experimental data for ^®°Sm [panel 
(a)] oscillate with J but become zero around J « S'*". 
Though larger deviations are observed in our calcula¬ 
tions [panel (b)] their global trend resembles the exper¬ 
imental one. Both theoretically and experimentally, the 
deviation from S{J) = 0 in ^^®Sm is more pronounced 
than for This suggests a deviation from the 

ideal alternating-parity band behavior, and also corre¬ 
lates well with the behavior of the Gogny-DIM MFPESs 
(see, Fig.[^. 


The reduced transition probabilities B{E3'^i — )• 0)'’) 
and B{E1;1^ —)• 0)*') are compared in Fig. ^ with the 
experimental data [TTllSOj . For both isotopic chains, the 
predicted E3 transition rates [panels (a) and (b)] exhibit 
a weak dependence on the neutron number with a max¬ 
imum at N = 88 — 90. The down-sloping tendency in 
the theoretical (IBM) E3 values observed in the heav¬ 
ier isotopes is consistent with the experiment though a 
smoother change with neutron number is found for Sm 
isotopes. On the other hand, the FI transition rates [pan¬ 
els (c) and (d)] increase with increasing neutron number 
which agrees quite well with the experiment, exception 
made of ^^®Sm. The overall trend also agrees well with 
the one found in previous IBM and GCM [19] calcu¬ 
lations. Note that the discrepancy of the IBM rates with 
the experimental ones are partly a consequence of the 
particular choice of the IBM effective charges. No effec¬ 
tive charges are needed within the GCM framework [19] 
as all the nucleons are considered in the wave functions. 

B. Spectroscopy of the nuclei and 

The low-lying spectrum of is compared in Fig.[^ 

with the available experimental excitation energies m 
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FIG. 5. (Color online) The same as in Fig. but for the lowest-lying odd-spin negative-parity states up to = 9 
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FIG. 6. (Golor online) The excitation energies of the 1)“ states predicted with the mapped IBM Hamiltonian are compared with 
the ones obtained in the framework of two-dimensional GGM calculations [19] for Sm [panel (a)] and Gd [panel (b)] nuclei. In 
both methods, the Gogny-DIM parametrization has been used. The experimental energy levels are also included in the figure. 


The band assignment has been made according to the Approximate alternating parity bands can be seen with 
dominant E2 transition sequence. The IBM energies the level ordering 7~, 8^, 9“, 10'*", ...etc. 
are generally more stretched than the experimental ones. 

A noticeable deviation with respect to the experimen- 
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J 


FIG. 7. (Color online) Signature splitting S{J) of nuclei as a function of spin J. For more details, see the main 

text. 



FIG. 8. (Color online) Theoretical and experimental transition probabilities > 0)^) and B(i?l; Ij —>0)'') for 

i4®-i5®Sm and 44®“i5®Gd. The experimental data are taken from [TTHSO] . 


tal data is obtained for the /?-vibrational band-head. In fact, the experimental excitation energy of this Oj state 
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FIG. 9. Comparison of the low-energy spectrum predicted within the IBM framework for the nucleus ^®°Sm with the available 
experimental excitation energies U7\. 


is as small as the one for the 4^ state. However, in the 
calculations it is almost twice higher, suggesting a too 
limited IBM model space. On the other hand, for the 
quasi-7 band, with the K'^ = 2+ built on the 2^ state, 
our calculations predict the staggering (3+, 4+), (5+, 6+), 
etc. This reflects the lack of triaxiality in the present 
study. The inclusion of mean-field triaxiality as well as 
the relevant terms in the mapped IBM Hamiltonian could 
be useful to better describe the structure of the quasi-y 
band [52]. Work along these lines is in progress and will 
be reported elsewhere. 

The E2 and El transition rates obtained for ^^°Sm are 
compared with the experimental ones [77] in Tables [I] and 
[H] respectively. Most of the predicted E2 values agree 
reasonably well with the experiment. Note that our cal¬ 
culations account for the K = Q~ band, built on the 3]" 
state, with strong E2 transitions. Nevertheless, large dis¬ 
crepancies are also found for some inter-band transitions. 
For example, the 0^ —>■ 2'^ strength is considerably un¬ 
derestimated. Stronger inter-band E2 transitions suggest 
a significant mixing between different intrinsic configura¬ 
tions. Indeed, a recent experiment has suggested a com¬ 
plex shape coexistence in ^®^Sm [5|. Within this context, 
an IBM model space larger than the one considered in the 
present study may be required. A configuration mixing 
associated with intruder states |53| could also be intro¬ 
duced to better describe a transitional nucleus like ^^°Sm. 
Another alternative could be the inclusion of triaxiality 
to better constrain the form of the IBM Hamiltonian. 
Furthermore, the value B{E2] 4^ —>• 2+) = 0.039 W.u. is 
too small as compared with the experimental one [42(20) 


W.u]. A possible reason may be that the Ot states as 
well as the ones built on it might not be well described 
by the present calculations. 

The calculated B{E1) values in Table |h] reveal rather 
strong transitions (starting around the J > 5“) from 
the states of odd-J negative-parity AT = 0]" to those 
of the even-( J — I) positive-parity ground-state bands. 
This fact, as well as the increasing B{E1;J^^^_ —>■ 
(J — value, as a function of J, signals the exis¬ 

tence of an alternating parity band in ^®'^Sm. Neverthe¬ 
less, we do not consider the B{E1) value obtained in the 
present calculation to be conclusive, mainly because of 
the lack of the p-boson effect in our framework. Indeed, 
as already pointed out in previous phenomenological m 
and microscopic |8T] studies on octupole-deformed nu¬ 
clei, the description of these El transitions in the IBM 
framework could be improved by explicitly including the 
p boson in the model space or by extending the form of 
the El operator so as to absorb the p-boson effect in the 
sdf space. 

The low-lying spectrum of ^®®Gd, shown in Fig. 
exhibits an overall agreement with the available experi¬ 
mental data for the lowest-lying positive- and negative- 
parity bands. The 1]” state of the lowest negative-parity 
band is assigned as the band-head of the — 0]" and 
the = ir bands in the present calculation and in the 
NNDC compilation [77], respectively. 

In our calculations, the two lowest-lying positive- 
parity, K'^ = 0]/" and 2+, bands are comprised of states 
with Uf ~ 0.02 and 0.06 ^ n/ < 0.09, respectively, sug- 
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FIG. 10. Same as in Fig. [^but for the nucleus ^®®Gd. 


gesting that they are almost pure positive-parity bands. 
On the other hand, the states in the band built on the 
Oj (0^) state are of two-/ boson (equivalently double 
octupole phonon) nature with {fif) « 2. The side-band 
energies, especially for those states in the positive-parity 
/I-vibrational and quasi-7 bands, are overestimated con¬ 
siderably, for similar reasons as in the ^®°Sm case. 

In Tables. Ill and IV we have compared some relevant 
E2 and El transition rates with the experimental ones 
[77]. Many of the calculated E2 transition rates agree 
well with the data. Again a noticeable deviation is ob¬ 
served for the 4^ —>■ 2+ transitions, probably for the same 
reason as in the ^®°Sm case. We note that the lifetime of 
the experimental state adapted in |77| has nearly 

25% of uncertainty, and that, for that reason, the error 
bars for the reduced E2 transitions 4 


K=i- 


^K=i- 


5y _.|_ —>■ shown in Ref. as well as in Ta¬ 


ble could be corrected. From Table [IV| one concludes 
that our model gives a reasonable description of the El 
transitions associated to states in the = 0/ band 
whose energies are described rather nicely as well (see 
Fig. 10). However, our model does not account for the 
El transitions associated to the K'^ = 1/ band. 


C. Excited 0”*^ states 

Experimentally, many excited 0“*" states have been 
identified in the low-energy excitation spectrum of ^^®Gd. 
The previous phenomenological calculation within the 
spdf IBM model m showed that such a large number 
of excited O’*" states at relatively low energy can be de¬ 


scribed if the octupole degrees of freedom is taken into 
account, and many of the 0^ states have been attributed 
to the coupling of two octupole phonons. Meanwhile, 
the emergence of a large number of low-energy excited 
O'*' states can be a good signature of a quantum phase 
transition [82] . 


To address the nature of the 0+ states resulting from 
the mapped sdf IBM Hamiltonian, we show in Fig.fTTjthe 
energy distribution (or level scheme) of the lowest fifteen 
0+ states and the corresponding average values of the 
/-boson number operator (h/) for the [from 

panel (a) to panel (f)] and i48-158Qj^ [from panel (g) 
to panel (1)] nuclei. In the [panel (a)] and ^^®Gd 

[panel (g)] cases those states with an energy higher than 8 
MeV are not shown. In all the nuclei, the 0+ ground-state 
is predominantly composed of positive-parity (s and d) 
bosons as (h/) < 0.5. In both isotopic chains, for many 
of the nuclei with N ^ 90, (hf) « 2 for the Oj state, 
suggesting its double-octupole phonon nature. Moreover, 
many other 0^ states are also formed by the coupling of 
positive- and negative-parity (octupole) bosons. For both 
Sm and Gd chains, the 0+ states become more populated 
in lower-energy region for the heavier isotopes, where the 
quadrupole-octupole coupling becomes more enhanced. 
Particularly in the Gd isotopes, the level scheme for the 
0+ states becomes most compressed around ^®^Gd [panel 
(i)] or ^^'^Gd [panel (j)], where the corresponding poten¬ 
tial energy surface is noticeably soft both in (32 and 
deformations [see. Figs. SiJ)]. 
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FIG. 11. (Color online) Energy distribution of the theoretical lowest fifteen 0’*' states and expectation value of the /-boson 
number operator (n/) for the considered Sm and Gd isotopes. Note that, concerning [panel (a)] and ^"^®Gd [panel (g)], 

the 0"*" states with an energy higher than 8 MeV are not shown. 


D. Correlation energy 


V. SUMMARY 


In this section, we discuss the correlation energies de¬ 
fined as [m 150] 

Ecorr = El,;-^-E{0t). (14) 

where i?HFB represents the HFB ground-state energy 
and i?(0[j'') the one for the 0[j^ state. The IBM corre¬ 
lation energies are depicted in Fig. together with the 
ones obtained in previous two-dimensional Gogny-DIM 
GCM calculations. Results are shown in panel (a) for 
i46-i56gj^ and in panel (b) for ^^®“^®®Gd. Though the 
correlation energies are different in both approaches, the 
largest values of Ecorr are obtained for the lighter nuclei 
with N ^ 88 which are rather soft in the P 2 and de¬ 
grees of freedom. This confirms that correlations beyond 
the mean-field approach can become significant in soft 
nuclear systems. 


In summary, we have carried out spectroscopic calcula¬ 
tions aimed to describe the quadrupole and octupole col¬ 
lective states in Sm and Gd isotopes. Our starting point 
was a set of Q20 — Q30 constrained HFB calculations, with 
the DIM parametrization of the Gogny effective interac¬ 
tion, used to produce a potential energy surface. This 
potential energy surface is then used to obtain the param¬ 
eters of an IBM Hamiltonian including s, d and / bosons. 
Spectral properties of both positive- and negative-parity 
states associated to the reflection symmetric and asym¬ 
metric shapes, respectively, are obtained after diagonal- 
ization of the IBM Hamiltonian. The parameters of the 
IBM Hamiltonian are determined by mapping the Gogny- 
HFB mean-field energy surface onto the corresponding 
energy expectation value of the boson condensate state. 

The systematics of the energy spectra and transi¬ 
tion rates, associated to both positive- and negative- 
parity yrast states, points to the onset of notable oc- 
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Q I _I_I_ ^ _I_I_ ^ _ I _ ^ _I_I_ ^ _I_ ^ _ 

84 86 88 90 92 94 84 86 88 90 92 94 

Neutron number Neutron number 

FIG. 12. (Color online) The correlation energies obtained from the IBM and the two-dimensional (2D) GCM for 
and 148-158 isotopes. 


tupole correlation around N « 88, characterized by the 


TABLE I. Theoretical and experimental ^7] B(E2) transi¬ 
tions for (in Weisskopf units). For details, see the main 

text. 


7^ 

Jf 

B{E2)theoT. 

B{E2)e^pt. 

2^ 

of 

79 

57.1(13) 

4^ 

2^ 

112 

110(17) 

6^ 


120 

1.5xl0+^(5) 

8^ 

et 

117 

1.7x10+2(9) 

oj 

2^ 

10 

53(5) 


0^ 

1.78 

0.8ll2? 


0^ 

24 

1.1x102 Ig 


2^ 

0.038 

- 


4^ 

3.92 

- 

2+ 

0^ 

1.18 

2.1(15) 


0^ 

10.2 

9.174) 



21 

- 



4.09 

- 


4^ 

0.064 

7(3) 

3^ 

2^ 

11.1 

- 


2? 

0.86 

- 


2+ 

^/3 

11.2 

1.9x10+2(9) 


2+ 

^7 

0.039 

42(20) 


3^ 

0.34 

- 

4^ 

^7 

2^ 

0.14 

1.4(7) 


2 + 

^13 

2.1 

4.1(21) 


2+ 

42 

- 

ir 

3 k = o - 

109 

- 


3^" n- 

70 

- 

7L ’ 

57- . 

84 


K = Q 

K=0 



/Sa-soft energy surfaces (Fig. [^, and the corresponding 
negative-parity band lowering in energy with respect to 
the positive-parity ground-state band (Fig. [^. From 
TV ^ 90 on, the potential energy surface no longer ex¬ 
hibits Ps softness, and the corresponding negative-parity 
band is pushed up in energy with respect to the ground- 
state band. The mean-field P 2 P 3 energy surface (Fig. 1|, 
the derived parameters in the sdf Hamiltonian (Fig. 31, 
the resultant energy levels (Figs. and and transi¬ 
tion rates (Fig. correlate very well with each other 
in systematics with the number of valence nucleons. In 
addition, the spectroscopic properties resulting from the 


TABLE 11. Same as in Table [I] but for the El transitions (in 
10-® W.U.). 


7^ 

'^i 

Jf 

B(E'l)theor. 

i3(i51)expt. 


Ot 

1.1 

L4^^ 

2 ^ 

0.13 

2.911^ 

3ir=o- 

2 ^ 

2.5 

cr+4 

^-3 

4^ 

1.8xl0"'‘ 

c+4 

^-3 

4j 

3^=0- 

2.6 


4| 


0.24 

- 

3if-o- 

0.087 

0.27(13) 



0.54 

0.9(5) 


4^ 

4.2 

- 

6 ^ 


0.027 

- 


6 ^ 

5.8 

- 

8^ 

'^K=0- 

0.15 

- 

10+ 

8^ 

7.3 

- 

10+ 

0.46 

- 


9.0 

- 
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model turn out to be generally in a reasonable agreement 
with the systematics of the available experimental data, 
and also to be consistent with the previous GCM calcu¬ 
lation (Figs. and starting from the common Gogny 

parametrization DIM uni¬ 
on the other hand, an in-depth analysis of the energy 
spectra and the E2 and El transition rates in the two 
characteristic cases, (32- and /Sa-soft nucleus ^^°Sm and 
strongly f }2 deformed nucleus ^®®Gd, has revealed that 
an improvement of the model is required so as to give a 
better description not only of the yrast states but also 
of the non-yrast states. For example, our model in its 
current version is not able to describe well the band- 
head of side bands, particularly that of the /3-vibrational 
[K^ = 0+) band (Figs. § and . A possible reason 
could be that the model space used for the present work 
might be rather limited to handle such a complex nu¬ 
clear structure. This would require the extension of our 


TABLE III. Same as in Table[i]but for the nucleus ^®®Gd. 


7^ 

Jl 


5(£^2)expt. 

2 ^ 

0 ^ 

170 

198(6) 

4/ 

2 / 

241 

289(5) 

6 / 

4/ 

260 

- 

8 / 

6 / 

264 

3.3x10+2(3) 

2 + 

0 / 

3.9 

3.4(3) 


2 / 

7.7 

6.0(7) 


4/ 

0.59 

0.27(4) 


0 / 

0.14 

0.31(4) 


2 / 

0.088 

1.39(15) 


2 + 

0.88 

12.8 


2 ^ 

93 

455 

^K-l- 

iK.r- 

105 

- 


1^.0- 

146 

> 1.6 X 10+2 


^K-l- 

41 

781(14) 



189 

369(6) 


‘^K=0- 

159 

- 



142 

2.09x10+2(3) 


model space to include configuration mixing specific to 
the intruder state and/or to introduce triaxial degrees of 
freedom. In addition, the model has failed in reproducing 
some of the El properties, especially for those associated 
to the states in non-yrast negative-parity band (Eig. 101. 
Several solutions have been proposed that could help to 
fix the problem: extension of the El operator to include 
higher-order terms; explicit inclusion of p boson in the 
model space. Improving the description of these proper¬ 
ties will be a topic of future study. Significance of the 
p-boson effect in the El excitation observed in rare-earth 
nuclei has been addressed in [53] , though in the different 
context of a clustering. 


We have also analyzed the wave function content of 
some lower-lying excited O’*" states for the considering 
nuclei, and found that in many of the nuclei considered, 
the 0^ states can be the consequence of the coupling of 
two-octupole phonons. This could be a possible explana¬ 
tion for the large number of low-energy excited O'*" states 
found in rare-earth nuclei. 


TABLE IV. Same as in Tablej^but for the El transitions (in 
units of 10“^ W.n.) in i®®Gd. 


Jj Jf B(El)theoT. I3(ifl)expt. 


Ik.i- 

0 ^ 

1.5 

0.098443(4) 


2 / 

2.8 

0.096515(6) 

^K-O- 

0 / 

4.3 

3.5(12) 


2 / 

2.3 

6.4(21) 

^K-1- 

2 / 

0.35 

0.33(10) 


4/ 

3.4 

0.29(8) 

CO 

1 

II 

o 

1 

2 / 

6.8 

> 1.1 


4/ 

0.74 

> 1.5 

2^.1- 

2 / 

4.5 

< 0.078 


4/ 

4.5 

0.090628(4) 


ACKNOWLEDGMENTS 


K. N. acknowledges the support by the Marie Curie 
Actions grant within the Seventh Eramework Program 
of the European Commission under Grant No. PIEF- 
GA-2012-327398. The work of LMR is supported in part 
by Spanish MINECO grants Nos. FPA2012-34694 and 
FIS2012-34479 and by the Consolider-Ingenio 2010 pro¬ 
gram MULTIDARK CSD2009-00064. 














16 


[11 P. A. Butler and W. Nazarewicz, Rev. Mod. Phys. 68, 
349 (1996). 

[2] P. A. Butler and W. Nazarewicz, Nucl. Phys. A 533, 249 
(1991). 

[3] L. P. Gaffney, P. A. Butler, M. Scheck, A. B. 
Hayes, F. Wenander, M. Albers, B. Bastin, C. Bauer, 

A. Blazhev, S. Bonig, et al.. Nature (London) 497 , 199 
(2013). 

[4] W. Urban, R. Lieder, W. Cast, G. Hebbing- 

haus, A. Krmer-Flecken, K. Blume, and H. Hbel, 
Physics Letters B 185 , 331 (1987), ISSN 0370- 

2693, URL http://www.sciencedirect.com/science/ 
article/pii/0370269387910094 

[5] W. Urban, R. Lieder, J. Bacelar, P. Singh, D. Alber, 
D. Balabanski, W. Gast, H. Grawe, G. Hebbinghaus, 
J. Jongman, et al., Physics Letters B 258, 293 (1991), 
ISSN 0370-2693, URL http://www.sciencedirect.com/ 
science/article/pii/037026939191088D 

[6] P. E. Garrett, W. D. Kulp, J. L. Wood, D. Bandyopad- 
hyay, S. Ghoudry, D. Dashdorj, S. R. Lesher, M. T. McEl- 
listrem, M. Mynk, J. N. Orce, et ah, Phys. Rev. Lett. 
103, 062501 (2009). 

[7] R. F. Gasten and N. V. Zamfir, Phys. Rev. Lett. 
87, 052503 (2001), URL http://link.aps.org/doi/10. 
1103/PhysRevLett.87.052503 

[8] F. lachello, Phys. Rev. Lett. 87, 052502 (2001), URL 
http://link.aps.org/doi/10.1103/PhysRevLett.87. 
052502 

[9] S. R. Lesher, A. Aprahamian, L. Trache, A. Oros- 
Peusquens, S. Deyliz, A. Gollwitzer, R. Hertenberger, 

B. D. Valnion, and G. Graw, Phys. Rev. C 66, 
051305 (2002), URL http://link.aps.org/doi/10. 
1103/PhysRevC.66.051305 

[10] N. V. Zamfir, J.-y. Zhang, and R. F. Gasten, Phys. Rev. 
G 66, 057303 (2002), URL http://link.aps.org/doi/ 
10.1103/PhysRevC.66.057303 

[11] S. Marcos, H. Flocard, and P. H. Heenen, Nucl. Phys. A 
410 , 125 (1983). 

[12] W. Nazarewicz, P. Glanders, 1. Ragnarsson, J. Dudek, 
G. A. Leander, P. Moller, and E. Ruchowsa, Nucl. Phys. 
A 429 , 269 (1984). 

[13] W. Nazarewicz and P. Glanders, Nucl. Phys. A 441 , 420 
(1985). 

[14] P. Bonche, P.-H. Heenen, H. Flocard, and D. Vautherin, 
Phys. Lett. B 175 , 387 (1986). 

[15] P. Bonche, in The Variation of Nuclear Shapes, edited by 
J. D. Garrett (World Scientific, Singapore, 1988), p. 302. 

[16] J. L. Egido and L. M. Robledo, Nucl. Phys. A 524, 65 
(1991). 

[17] L. M. Robledo, M. Baldo, P. Schuck, and X. Viiias, Phys. 
Rev. C 81, 034315 (2010). 

[18] L. M. Robledo and G. F. Bertsch, Phys. Rev. G 84 , 
054302 (2011). 

[19] R. Rodrfguez-Guzman, L. M. Robledo, and P. Sarriguren, 
Phys. Rev. G 86, 034336 (2012). 

[20] L. M. Robledo and P. A. Butler, Phys. Rev. C 88, 051302 
(2013). 

[21] O. Scholten, F. lachello, and A. Arima, Ann. Phys. (NY) 
115 , 325 (1978). 

[22] J. Engel and F. lachello, Nucl. Phys. A 472, 61 (1987). 

[23] T. Gtsuka and M. Sugita, Phys. Lett. B 209, 140 (1988). 


[24] D. Kusnezov and F. lachello, Phys. Lett. B 209, 420 
(1988). 

[25] P. D. Gottle and N. V. Zamfir, Phys. Rev. C 58, 1500 
(1998). 

[26] P. G. Bizzeti and A. M. Bizzeti-Sona, Phys. Rev. G 70 , 
064319 (2004). 

[27] D. Bonatsos, D. Lenis, N. Minkov, D. Petrellis, and 
P. Yotov, Phys. Rev. G 71 , 064309 (2005). 

[28] D. Lenis and D. Bonatsos, Phys. Lett. B 633, 474 (2006). 

[29] P. G. Bizzeti and A. M. Bizzeti-Sona, Phys. Rev. G 77 , 
024320 (2008). 

[30] P. G. Bizzeti and A. M. Bizzeti-Sona, Phys. Rev. G 81, 
034320 (2010). 

[31] R. V. Jolos, P. von Brentano, and J. Jolie, Phys. Rev. C 
86, 024319 (2012). 

[32] N. Minkov, S. Drenska, M. Strecker, W. Scheid, and 
H. Lenske, Phys. Rev. C 85, 034306 (2012). 

[33] P. G. Bizzeti and A. M. Bizzeti-Sona, Phys. Rev. G 88, 
011305 (2013). 

[34] F. lachello and A. D. Jackson, Phys. Lett. B 108, 151 
(1982). 

[35] H. J. Daley and F. lachello, Ann. Phys. (NY) 167, 73 
(1986). 

[36] T. M. Shneidman, G. G. Adamian, N. V. Antonenko, 
R. V. Jolos, and W. Scheid, Phys. Lett. B 526, 322 
( 2002 ). 

[37] W. Zhang, Z. P. Li, S. Q. Zhang, and J. Meng, Phys. 
Rev. C 81, 034302 (2010). 

[38] B.-N. Lu, E.-G. Zhao, and S.-G. Zhou, Phys. Rev. G 85, 
011301 (2012). 

[39] B.-N. Lu, J. Zhao, E.-G. Zhao, and S.-G. Zhou, Phys. 
Rev. C 89, 014323 (2014). 

[40] M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. 
Phys. 75, 121 (2003). 

[41] T. H. R. Skyrme, Nucl. Phys. 9, 615 (1958). 

[42] D. Vautherin and D. M. Brink, Phys. Rev. C 5, 626 
(1972). 

[43] J. Decharge and M. Girod and D. Gogny, Phys. Lett. B 
55, 361 (1975). 

[44] D. Vretenar, A. V. Afanasjev, G. Lalazissis, and P. Ring, 
Phys. Rep. 409 , 101 (2005). 

[45] T. Niksic, D. Vretenar, and P. Ring, Prog. Part. Nucl. 
Phys. 66, 519 (2011). 

[46] R. Rodrfguez-Guzman, J. L. Egido, and L. M. Rob¬ 
ledo, Nuclear Physics A 709 , 201 (2002), ISSN 0375- 
9474, URL http://www.sciencedirect.com/science/ 
article/pii/S0375947402010199 

[47] P. Ring and P. Schuck, The nuclear many-hody problem 
(Berlin: Springer-Verlag, 1980). 

[48] L. M. Robledo and R. R. Rodrguez-Guzman, Journal 
of Physics G: Nuclear and Particle Physics 39, 105103 
(2012), URL http: //stacks . iop. org/0954-3899/39/i= 
10/a=105103 

[49] K. Nomura, N. Shimizu, and T. Otsuka, Phys. Rev. Lett. 
101, 142501 (2008). 

[50] K. Nomura, N. Shimizu, and T. Otsuka, Phys. Rev. C 
81, 044307 (2010). 

[51] K. Nomura, T. Otsuka, N. Shimizu, and L. Guo, Phys. 
Rev. G 83, 041302 (2011). 

[52] K. Nomura, N. Shimizu, D. Vretenar, T. Niksic, and 
T. Otsuka, Phys. Rev. Lett. 108, 132501 (2012). 



17 


[53] K. Nomura, R. Rodriguez-Guzman, L. M. Robledo, and 

N. Shimizu, Phys. Rev. C 86, 034322 (2012). 

[54] K. Nomura, D. Vretenar, and B.-N. Lu, Phys. Rev. C 88, 
021303 (2013). 

[55] K. Nomura, D. Vretenar, T. Niksic, and B.-N. Lu, Phys. 
Rev. C 89, 024312 (2014), URL http://link.aps.org/ 
doi/10.1103/PhysRevC.89.024312 

[56] S. Goriely, S. Hilaire, M. Girod, and S. Peru, Phys. Rev. 
Lett. 102 , 242501 (2009). 

[57] R. Rodriguez-Guzman, P. Sarriguren, L. M. Robledo, 
and S. Perez-Martin, Phys. Lett. B 691 , 202 (2010). 

[58] R. Rodriguez-Guzman, P. Sarriguren, and L. M. Robledo, 
Phys. Rev. G 82, 044318 (2010). 

[59] R. Rodriguez-Guzman, P. Sarriguren, and L. M. Robledo, 
Phys. Rev. G 82 , 061302 (2010). 

[60] S. A. Giuliani, L. M. Robledo, and R. Rodriguez- 
Guzman, Phys. Rev. G 90, 054311 (2014), URL http: 
//link.aps.org/doi/10.1103/PhysRevC.90.054311 

[61] J. F. Berger, M. Girod, and D. Gogny, Nucl. Phys. A 
428 , 23 (1984). 

[62] Notel, we equally use the multipole moment values Qio 
and deformation parameters Pi to talk about deforma¬ 
tion. 

[63] T. Otsuka, A. Arima, and F. lachello, Nucl. Phys. A 309 , 
1 (1978). 

[64] M. Sugita, T. Otsuka, and P. von Brentano, Phys. Lett. 
B 389 , 642 (1996). 

[65] A. F. Barfield, B. R. Barrett, J. L. Wood, and 

O. Scholten, Ann. Phys. 182 , 344 (1988). 

[66] J. N. Ginocchio and M. W. Kirson, Nucl. Phys. A 350, 
31 (1980). 

[67] N. V. Zamfir and D. Kusnezov, Phys. Rev. C 63, 054306 

( 2001 ). 

[68] M. Babilon, N. V. Zamfir, D. Kusnezov, E. A. Mc- 
Cutchan, and A. Zilges, Phys. Rev. G 72 , 064302 (2005). 


[69] D. J. Thouless and J. G. Valatin, Nucl. Phys. 31, 211 
(1962). 

[70] H. Schaaser and D. M. Brink, Nucl. Phys. A 452, 1 
(1986). 

[71] D. Kusnezov, computer program OCTUPOLE (unpub¬ 
lished) . 

[72] W. Long, J. Meng, N. V. Giai, and S.-G. Zhou, Phys. 
Rev. C 69, 034319 (2004). 

]73] T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. G 78, 
034318 (2008). 

[74] T. Otsuka, Phys. Rev. Lett. 46, 710 (1981), URL http: 
//link.aps.org/doi/10.1103/PhysRevLett.46.710 

[75] T. Otsuka and J. N. Ginocchio, Phys. Rev. Lett. 55, 
276 (1985), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.55.276 

[76] F. lachello and A. Arima, The interacting boson model 
(Gambridge University Press, Cambridge, 1987). 

[77] Brookhaven National Nuclear Data Center, 

http://www.nndc.bnl.gov. 

[78] F. R. Metzger, Phys. Rev. C 14, 543 (1976). 

[79] T. Kibedi and R. H. Spear, At. Data and Nucl. Data 
Tables 80, 35 (2002). 

[80] H. Pitz, U. Berg, R. Heil, U. Kneissl, R. Stock, 

C. Wesselborg, and P. V. Brentano, Nuclear 

Physics A 492, 411 (1989), ISSN 0375-9474, URL 

http://www.sciencedirect.com/science/article/ 
pii/0375947489904053 

[81] T. Otsuka, Phys. Lett. B 182, 256 (1986). 

[82] D. Meyer, V. Wood, R. Casten, C. Fitzpatrick, G. Graw, 

D. Bucurescu, J. Jolie, P. von Brentano, R. Hertenberger, 
H.-F. Wirth, et ah. Physics Letters B 638, 44 (2006), 
ISSN 0370-2693, URL http://www.sciencedirect.com/ 
science/article/pii/S0370269306005594 

[83] M. Spieker, S. Pascu, A. Zilges, and F. lachello, Phys. 
Rev. Lett. 114, 192504 (2015), URL http://link.aps. 
org/doi/10.1103/PhysRevLett.114.192504 


